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We show that an inflationary slow-roll potential can be derived as an IR limit of the non- 
perturbative exact renormalisation group equation for a scalar field within the mean-field approxi¬ 
mation. The result follows without having to specify a Lagrangian in the UV, which we take to be 
somewhere below the Planck scale to avoid discussing quantum gravity effects. We assume that the 
theory contains a scalar mode with suppressed coupling to other fields, and that higher derivative 
couplings are suppressed. In this framework the exact RG equation becomes a one-dimensional 
Schrodinger equation, which we solve. The effective IR potential is then dominated by the eigen¬ 
states of the RG Hamiltonian with the highest eigenvalues. We find that these potentials can 
generically give rise to slow-roll inflation, which is fully consistent with recent observations. As 
an example of how the proposed renormalisation group procedure works, we perform an explicit 
calculation in the <f ) 4 theory in an appendix. 
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I. INTRODUCTION 

Inflation [1-3] is an exponentially fast expansion in 
the early universe that is claimed to resolve a number 
of problems of the standard Big Bang cosmology (e.g. 
the horizon, flatness, smoothness, relic problems). It 
also provides a mechanism for generation of the den¬ 
sity perturbations [4] that have left an imprint on the 
cosmic microwave background and grown into the ob¬ 
served large-scale structure of the universe [5]. Models 
of inflation are typically discussed within the slow-roll 
paradigm, where an inflaton field evolves along a nearly 
flat potential to its minimum, during which it sources a 
quasi-exponential growth of the cosmological scale fac¬ 
tor. Hence, in the slow-roll regime, the spatial and time 
derivatives are negligible compared to the approximately 
constant vacuum energy. 

Many potentials driving the inflationary fields have 
been proposed that are either phenomenological or in¬ 
spired by various UV theories [6]. In this paper, we 
propose a different approach to deriving a scalar field 
theory suitable for the slow-roll inflation: We will ar¬ 
gue that inflation can be understood as a very generic, 
non-perturbative prediction of the exact renormalisation 
group (ERG), in a way that is largely insensitive to the 
details of the UV physics which we take to be at a scale 
somewhere below the Planck scale, Mpy < Mpi , in order 
to ignore quantum gravity effects. We add that the obser¬ 
vation that renormalisation group dynamics can produce 
experimentally viable inflation scenarios has been noted 
in the literature before [7, 8]. Moreover, exact RG has 
also been applied in the context of inflation before, see 
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e.g. [9-11]. 

To have a valid potential for the description of infla¬ 
tion, we must ensure that the effective theory is valid 
below some scale M, where M is roughly on the order of 
the inverse Hubble radius at the start of inflation. This 
is the IR of our theory. Thus, 1/M 4 will be approxi¬ 
mately the volume of the patch of spacetime of the re¬ 
sulting effective theory, which we assume resides within 
the Hubble radius. We will assume that the UV theory 
at Muv » M contains a scalar mode (j) that couples 
very weakly to other fields. Hence, during the RG flow 
from Muv to M, the corrections to the effective poten¬ 
tial of such a scalar from the couplings to other UV fields 
would remain small. In the RG analysis, <p can therefore 
be treated independently. 

The second assumption will be the validity of the 
mean-held approximation (MFA) near the scale M. 1 
That is, we will treat the constant mode of the theory 
separately when we do our IR analysis. In the paper, 
we will give an argument for the plausibility of this as¬ 
sumption, given that a held can be expanded around its 
constant value. We work in the local potential approx¬ 
imation, where higher derivative couplings are ignored. 
The argument will then rely on the particular form of the 
ERG equation. As a result, we will only be interested in 
the constant mode of the theory, which will be sufficient 
to describe the effective potential. Without justihcation, 
this assumption may at hrst seem peculiar, but we note 
that for inflation to begin, the inhaton held has to be 
sufficiently smooth to overcome the gradient energy pre¬ 
venting the exponential expansion. Indeed, the MFA is 
a common feature of many other inhationary scenarios, 
where the held is assumed constant at the scale where 
inhation begins [15]. 

Lastly, we will also ignore effects of gravity in the ERG. 
We start at a scale Muv somewhere below the Planck 
scale, so that quantum gravity effects can be ignored. For 
the most part of the RG how, we are at an energy scale 
where the vacuum energy is negligible in comparison, and 
we can work in Minkowski space where gravitational ef¬ 
fects are unimportant. Indeed, the de Sitter value of 
the vacuum energy of the universe is roughly comparable 
with the scale at which inhation begins, which is where 
we end our how. Flows beyond this point would require 
us to take gravitational effects into account [12-14], 

After introducing the ERG equation relevant for our 
setup and justifying the validity of the MFA, we will show 
that the MFA of the ERG equation gives rise to a one¬ 
dimensional Schrddinger-type equation for the resulting 
effective theory, which can be solved exactly. This is sim¬ 
ilar to the stochastic approach of [16, 17] as well as the 
study of renormalisation group in quantum mechanics 
[18]. In essence, we will solve a simple quantum mechan¬ 
ical problem in which the resulting “wave-function” will 


1 For other recent implementations of the ERG within the MFA, 
in particular in de Sitter space, see [12-14]. 


correspond to the effective theory, or potential, for the 
constant mode. Because in QFT the theories need not be 
normalisable functions of the field variable, we will only 
insist that the resulting potential is bounded from below. 
What is meant by normalisable and un-normalisable so¬ 
lutions will become clear below. Our analysis will then 
result in non-perturbative IR effective potentials with 
shapes suitable for inflation. It will give predictions that 
are fully consistent with current observations. We note 
that spectral methods have been applied in the context 
of ERG before, see e.g. [19-23] and references therein for 
more details. 

At the end of the paper, we provide two appendices: 
A is devoted to an ERG analysis of scalar field theories, 
in particular the <?!> 4 -theory, which is an example of a nor¬ 
malisable theory. B is devoted to solving the Schrodinger 
equation in momentum space, where the solution takes a 
rather simple form. 


II. EXACT RENORMALISATION GROUP 
ANALYSIS 


A. The ERG equation 


Let us begin by considering a UV quantum held theory 
of the early universe that contains some set of fields at 
the UV scale Muv < Mpi, including a scalar mode <f>. 
By 4>, we collectively denotes the remaining UV fields. 
As the theory runs between Muv and the IR scale M, 
we assume that </> and $ interact very weakly so that the 
IR theory for p < M takes the form, 


Z= [ V(j)e- S * w [ £>$ e _s * [<61 . 
J M J M 


(1) 


We can thus focus only on the decoupled renormalisation 
group flow for the scalar <j> by following the procedure of 
the exact Wilsonian renormalisation group [24]. 2 The 
ERG equation for the decoupled scalar mode takes the 
form 



6 S^ SS 4 
54>p 


, SS, 


( 2 ) 


where a{t) depends on the choice of the cut-off and t ~ 
1/A denotes the RG time - an inverse of the cut-off scale 
A. 3 4 From an inflationary standpoint it is also convenient 
to measure the dimensionless field (f> in units of Mpi, 
rather than Muv - 4 


2 For reviews on the exact renormalisation group, see [25—27]. 

3 Since we are only interested in solving the theory in the IR, we 
do not include the rescaling of the theory back to the original 
UV scale — blocking from Mjjy to M suffices. 

4 Throughout this paper, we use the reduced Planck mass, A'lpi = 
2.44 x 10 18 GeV. 
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The ERG equation (2) is a non-linear equation. How¬ 
ever, if we instead consider the functional 

= e- a *M, (3) 

then Eq. (2) becomes a linear Schrodinger-type equation, 
dt'Z [t, <t>] = H'i [t, 4 >], (4) 

with a Hamiltonian operator in position space, 



This is the usual form of the Wilsonian ERG equation, 
written in terms of 4/. Similar equations have been de¬ 
rived in the literature, see e.g. [28, 29]. Note that one 
common feature of these equations is that the Hamilto¬ 
nian has a kinetic part S 2 /5(j ) 2 , together with a divergence 
term <f>S/S<j). 


B. Validity of the MFA 

Before performing a detailed analysis of Eq. (4), we 
pause to discuss more precisely why the MFA is a valid 
approximation of the ERG equation when we are only 
interested in the potential of the theory for an approx¬ 
imately constant field. Working in the local potential 
approximation, we find that the kinetic modes decouple 
faster as we move towards the IR, and we integrate them 
out from the theory. This leaves us with an effective po¬ 
tential for the constant mode. We should also note that 
our conventions are such that all fields and coordinates 
will be dimensionless throughout the RG analysis. 

We begin by expanding the functional T as 

4/[0]=4>°[«/> o ] + [ *p^[0o] 

Jp 

+ 9 / e t > P < l ) q^p,q[ ( l > o] + • ■ • > (6) 

z J{p,q} 

where 4>J = = 0, and 4/^ can be assumed to be 

symmetric in p and q. We ignore terms that are cubic or 
higher-order in momenta. This is closely related to the 
often used local potential approximation. From hereon, 
we will denote the constant mode as <fio = x. Moreover, 


V°[x] = e~ s °^\ (7) 

which gives the potential of the theory. Inserting this into 
Eq. (4), we can derive the following set of equations, 

+ f (a(t) + 2 q 2 ) ^l_ q , (8) 

Jq 

dt^), = + (a(t) + 2 p 2 ) (9) 

d t ^l q = Ho*l >q + 2 (a(t) +p 2 + q 2 ) (10) 


where 

Ho = a(t) (dl + xd x ) . 

Eqs. (8) and (10) can then be put in the form 


d t 


V ]/ 0 

vj / 2 

p,q 


Ho tr 

0 Hr,, 


\]/ 0 

\]> 2 

p,q 


where 

Hp.q = Ho + 2 (cu(t) + p 2 + q 2 ) , 
while the trace tr is given by 


( 11 ) 


( 12 ) 


(13) 


t" r (^,g) = [ (“(*) + 2 « 2 ) K~q- ( 14 ) 

Jq 

As t —> oo, we can assume that a(t) —> a becomes con¬ 
stant [27]. The solution of (12), as t —> oo, is therefore 
given by the highest eigenvalue eigenmodes of the matrix 


II = 




(15) 


This matrix is upper-triangular, so the eigenmodes fac¬ 
tor into the eigenmodes of Ho and H v . q . In particular, 
the dominant term in the potential will be given by the 
highest eigenmode of Ho- 

For completeness, let us also consider what happens to 
the kinetic terms in the expansion (6) as t oo. Note 
first that the term linear in (f > p , in Eq. (6), integrates to 
zero, and we will hence ignore it. Note also that in order 
for eigenmodes of Ho not to blow up at infinity, we will 
see below that the eigenvalue is required to be less than 
or equal to zero. Finally, note form (13) that the kinetic 
terms 4/^ will blow up quicker than 'F°, an d w iH hence 
integrate out sooner. This is true even for momenta close 
to zero due to the non-vanishing a. Ignoring higher-order 
kinetic terms, the action reads 


S[4>] = - log(tf°[a;]) - \ [ [*] + ■••, (16) 

z J {p,q} 

where & 2 q [x] = ^> 2 q [a;] /\E ,0 [:r]. Note that in order to 
have a positive-definite kinetic term, we require 'b 2 q [<^o] 
to be negative. If we also assume that saturates at 
the highest allowed mode, which is constant, i.e. 

K,q M = ~Cp, q , (17) 

then we can integrate out the momentum modes to get 
the effective action 

S{x) = S 0 (x) + ( S 0 (x ) + log(Gq ig )). (18) 

We thus see that the correction only ends up multiply¬ 
ing the potential with an overall constant and adding an 
irrelevant constant to the potential. Beyond that, the 
potential remains the same. Hence, we see that in the lo¬ 
cal potential approximation, we are justified in using the 
MFA and we will therefore ignore propagating modes for 
the remainder of the paper. 
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C. Solution of the RG flow 

We can now use the MFA to solve the ERG equation 
(4) for the patch of the universe of the size of 1/M 4 , i.e. 
in the IR regime of the theory. In the UV, we wish to 
remain as general as possible, so we do not specify the 
details of the theory. Let us then write the Euclidean 
partition function ( 1 ) with its initial theory specified in 
the UV as 

Z = j V^ uv [t,4>\. (19) 

In the MFA, an eigenmode of the Hamiltonian % evolves 
under the ERG equation (4) as 

*F(t, x) = e Et ^/(x), (20) 

where E are the eigenvalues of the RG time-independent 
equation 

W& = E'b. (21) 

Note that we are using = 'F(x). In analogy with quan¬ 
tum mechanics, we assume that 4 , [/y(t, x) can be ex¬ 
panded as 


In terms of the familiar quantum mechanical notation, 
we can use x = cp/Mpi for our initial UV field qb, where x 
is measured in units of Mpi as is convenient for compu¬ 
tations involving inflation. We also let p = —id x , which 
gives the Hamiltonian operator H: 

H=-(p 2 +p + \), ( 25 ) 

where we have set a = 1 without loss of generality. It 
is important to note that the solutions to (21) and (24), 
i.e. 4q and respectively, have the same eigenvalues, 
E^. In the MFA, the two Hamiltonians are related by 

H = e~ x 2 / 4 He x2 / 4 = d 2 + xd x . (26) 


The form of H in Eq. (25) is highly reminiscent of the 
quantum harmonic oscillator in imaginary (Euclidean) 
time. Indeed, by defining the ladder operators 


1 

d= -x + ip, 

(27) 

1 A 

a = -x - ip, 

(28) 

the Hamiltonian takes the form 

= -( 0*0 + 1 ). 

(29) 


'F uv{t,x ) = ( 22 ) 

i 

where Hq have eigenvalues Ei. In the IR, where t —► oo, 
the dominant contribution will therefore come from the 
highest eigenvalue solution, T, with max [Ei], that has a 
non-trivial overlap with the UV theory. The contribution 
of other eigenvalue solutions will decay exponentially fast 
with RG time so we expect the highest eigenvalue solu¬ 
tion to dominate even with little RG flow. Note that in 
Eq. (22), a functional integral over \Eq should be used 
instead of the sum when Ei are a continuous set. 

In Appendix A, we show an example of how the decom¬ 
position in Eq. (22) can be performed in the mean-field 
approximation for a UV theory with a “normalisable”, 
finite (path) integral of ’F uv(x ) over x. We use the sim¬ 
plest example of an interacting field theory: the </) 4 the¬ 
ory. In general, however, such an explicit computation 
of 7 i may be extremely challenging, especially for poten¬ 
tials which give formally divergent integrals. Hereon, we 
will therefore only assume that such a decomposition is 
possible and that we can treat the most general solution 
for the “ground state” theory (one with max/Ej]) as the 
dominant theory in the IR. See Appendix A for defini¬ 
tions and further details. 

It is convenient to introduce a new functional i/j so that 

4- = e - £ 2 / V- (23) 

The ERG eigenvalue equation in the MFA then becomes 
Hill = Eil). (24) 


Eq. (29) , which describes the renormalisation group evo¬ 
lution of the effective action thus only differs from the 
usual harmonic oscillator by an overall minus sign and 
the additive factor of 1 in place of 1 / 2 . 

For modes of i/j which tend to zero as x —> ±oo, the 
operator a* is the adjoint of a, a* = a '. To see this, let a 
and /? be modes that tend to zero at infinity and consider 

(aa,/3) = Jdx(aa)l3* = jdxa(a*/3)* = (a,a^/3), (30) 

where we have performed an integration by parts. H is 
negative definite for such modes, and the highest eigen¬ 
value is the “vacuum energy”, 

E 0 = E= -1, (31) 

which corresponds to the vacuum of the theory, ipo, 

ai/jo = 0. (32) 

The corresponding re-scaled “original” theory reads 

>Fo(a:) = C 0 e~ x2/4 i> 0 (x) 

= C 0 e~ x2 / 2 , (33) 

which precisely corresponds to the potential of a free the¬ 
ory. Besides the free theory, there exists a mode, ipo, with 
the same eigenvalue (31), 

4>o(x) = D 0 e~ x ~ /A e rfi ^ , (34) 

), (35) 


'Fo(z) = -Doe 1 / 4 ^o(a:) = D 0 e x /2 erfi (^= 
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were erfi(a;) denotes the imaginary error function. Note 
that while ip o(x) diverges asi-> ±oo, as expected for a 
un-normalisable mode of the harmonic oscillator, the re¬ 
scaled wave-function, 4/o, exhibits the correct behaviour 
and tends to zero at infinity. 4' 0 leads to a bounded 
potential, as 

S+M = “log ^o = ^, (36) 


where 


where Cq and Dq are arbitrary constants. Importantly, 
the IR effective theory is continuously connected to the 
free Gaussian IR limit, as is often the case in perturbative 
RG. 

If we lift the restriction of including the free Gaussian 
potential in our IR solution, then it becomes more nat¬ 
ural to consider the second, continuous class of theories 
for which the Eq. (24) (or Eq. (21)) gives 


4/(cc) =e 2 x2 


Go l-Fi 


(l+E' 1 x 2 \ 

V 2 ; 2 ; Tj 


Vo(<j>) = M 4 


log (D 0 ) + \x 2 


log 



• (37) 


+ D 0 x 



/2 + D 

3 a; 2 \l 

V 2 

’ 2’ 2 ) 


(41) 


In order to avoid considering theories of arbitrarily 
high eigenvalues, one reasonable assumption is that the 
UV theory has no overlap with modes 4b of which the 
eigenvalues would be higher than the “vacuum” energy 
D 0 , cf. Eq. (31). In physical terms, this means that we 
assume that the IR limit of 4'j/y is connected to the free 
theory - the state 'Jo¬ 
in the absence of such a restriction, it is clear that the 
eigenstates of (24), and thus also 'J, could have arbitrar¬ 
ily high eigenvalues. However, we want our dominant 
IR theory, denoted temporarily by 4' = 4/s, with some 
eigenvalue E, to tend to zero as x —> ±oo, as required by 
a stability condition that potentials need to be bounded 
from below. The eigenvalue equation (21) then reads 

1e = D4 ' e 

= d^E + xdx'&E • (38) 

If we multiply this equation by 4 'e and integrate over x, 
we get after integration by part, 

E^e = - f (dx'S’E ) 2 f (39) 

The main point is that integration by parts can only be 
performed for 4 ^e that tend to zero at infinity. Since the 
integrands on the right-hand-side of (39) are negative- 
definite, it follows that theories with bounded potentials 
must have E < 0. 

In the remainder of this work, we will focus on two 
classes of stable effective potentials: the IR effective the¬ 
ories with discrete and continuous spectra of Ei < 0. In 
the deep IR limit where t —> oo, the existence of a dis¬ 
crete spectrum with integral eigenvalues implies that the 
state with the highest eigenvalue, E 0 = —1, will domi¬ 
nate the first class of theories. In the continuous case, 
theories with 2£* close to 0 will dominate. 

We can now restate the above results in the following 
way: By assuming that 4'j/y has a non-zero overlap only 
with the states of integral eigenvalues E = —n for n £ N, 
the IR limit of the theory takes the generic form 

(40) 


T 0 (a:) = e ^ 


C 0 + D 0 erfi 


(y /2 


where iF\ (a; b; x 2 /2) is the confluent hypergeometric 
function of the first kind. The solution in Eq. (41) is 
valid for all E and reduces to (40) at E = — 1. 

A simple complete set of orthogonal polynomials (the 
Hermite polynomials) can be formed from the 4/ (x) func¬ 
tions in (41) by restricting the eigenvalues E to be nega¬ 
tive integers, E £ {—1, —2, —3,...}. In this way, we can 
form an infinite series of finite and “normalisable” contri¬ 
butions to the expansion of 'S’uv hr (22). Here, we define 
normalisable to mean that the integral over x of 4* (a;) is 
finite for some choice of E: 


dx 4'£(:r) = finite. 


(42) 


As already noted above, we explore this possibility in de¬ 
tail in Appendix A where we use this basis to decompose 
the ERG flow of normalisable scalar theories, in partic¬ 
ular the </i 4 theory. However, in general, we expect that 
the integrals over 4r(a:) need not converge to give us a 
physically acceptable effective potentials. In the remain¬ 
der of this paper, we will study such “un-normalisable” 
theories, for which the integrals diverge: 


dx ^e(x) —> oo. 


(43) 


In particular, the potentials that we will use as candi¬ 
dates for inflation will be of the un-normalisable type. 

It is important to keep in mind that both MFA so¬ 
lutions, (40) and (41), can be non-perturbative in the 
coupling constants of the original UV theory. However, 
the dependence of the effective potential on the couplings 
could only be computed if we had specified the UV the¬ 
ory and somehow computed the relevant overlap integrals 
that would reveal the true weights of 4/^ in 4>. We will not 
pursue this direction in this work and will only consider 
IR effective potentials with unspecified coupling constant 
dependence. 

Before we move on to considering the phenomenolog¬ 
ical implication of the two solutions, we note that our 
approach resembles that of the stochastic approach dis¬ 
cussed in [16, 17]. It is also similar to the approach of 
Halpern and Huang [30] and Periwal [31], where a lineari¬ 
sation of the ERG equation (2) was performed. Those 
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works showed that there are non-polynomial deforma¬ 
tions of the free Gaussian theory for which the Gaussian 
is IR-unstable. Such cases are also included in our anal¬ 
ysis, within the MFA. Furthermore, this gives credence 
to the inclusion of the extra mode ipo in addition to the 
Gaussian when considering theories of eigenvalue —1. It 
also prompts us to consider in detail theories with eigen¬ 
value greater than —1, which can be argued to be prefer¬ 
able over the free theory from the exact renormalisation 
group point of view. 


Finally, we also note that our non-perturbative re¬ 
sult differs from the conventional perturbative approach 
where the effective Lagrangian is expanded in the pow¬ 
ers of the field suppressed by some mass scale. Our ap¬ 
proach is also different to what is known as the Effective 
Field Theory of Inflation [32], where the background dy¬ 
namics sourced by a potential is assumed as given and 
the effective theory refers to the fluctuations around that 
background. 


III. INFLATION 
A. Inflationary potential 

Having found solutions that are likely to dominate the 
non-perturbative IR regimes of scalar theories, we now 
turn our attention to studying the phenomenological im¬ 
plications of such theories. In particular, we will show 
that solutions of Eq. (4) naturally lead to potentials ca¬ 
pable of sustaining inflation in the early universe. The 
physical requirement in the patch of spacetime where in¬ 
flation starts is that the temporal derivative of the field 
is small and the field is homogeneous [15, 33]. To de¬ 
scribe inflation, we are therefore primarily interested in 
the potentials for the constant mode of the theory, which 
given the argument in section IIB means that we can use 
our MFA results, where we effectively ignored the higher 
momentum modes in order to derive equation (24) and 
its solutions. 

We begin by restating the un-normalisable solution 
with a continuous spectrum of Ei , i.e. Eq. (41), as an 
effective potential of </>. Restoring the mass dimensions, 


YM 

M 4 



log 


i-Fi 


(l + E 1 1 
^ 2 ’ 2’2 



where Mpi is the reduced Planck mass.. The potential 
has an overall factor expressed in terms of some mass 
scale M. Furthermore, we have two integration constants 
c and C. Our universe has a very small cosmological 
constant, therefore the constant C needs to be fixed so 
that the vacuum energy at the minimum of the potential 
is zero, V m i n = 0. This constrains C to be a function of 
c, i.e. C(c). 

At the minimum of the potential, the mass of the scalar 
field is independent of the shape parameter c, 

EM^ 

®<tM(fimin) = > ■ ( 45 ) 

m Pl 

When c > 1, the potential and its derivatives do not 
depend on c (see Fig. 1). Hence, the potential has a 
stable shape and the exact value of c for c > 0(1) does 
not matter. 

We are particularly interested in the cases that dom¬ 
inate the IR regime of the RG flow, where E —> CP. 
For small \E\, the plateau region of the potential is not 
only flat, but also small (see Fig. 1). Such a potential 
can support slow-roll inflation as well as result in small 
amplitude of the scalar perturbations required by the ob¬ 
servations, thus alleviating the fine-tuning of M (see next 
section). 



(2 + E 3 1 
y 2 ’2’2 



■ C(c) 

M 4 ’ 
(44) 


Let us explore the plateau region. We can use the 
asymptotic expansion of the hypergeometric functions to 
obtain the leading behaviour that is logarithmic: 


V{4>) » M' 4 


1 + a log 



(46) 


This potential shape is reminiscent of the ‘Loop Inflation’ 
model 5 where the logarithmic dependence arises from the 
loop corrections that “spoil” the flatness of the inflation¬ 
ary potential (i.e. the jy-problem). This has been studied 
in the context of the F- and D-term inflation. 6 


B. The slow-roll analysis in theories with a 
continuous spectrum of E 

The background dynamics of the homogeneous scalar 
field in the FRW universe is governed by the continuity 


5 Frequently termed ‘Spontaneously broken SUSY’ model. 

6 For a review of these approaches, see Ref. [6] and references 
therein. 















0.004 
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and Friedmann equations: 

4> + 3H<j> + d^V = 0, (47) 

fl2 = 3 wX^ 2+vw )- m 

Here, H = a/a is the Hubble parameter and a(t) is the 
scale factor that depends on time, not the RG time that 
we used in Sec. II. It follows then that the accelerated 
expansion of the universe (a/a > 0) is achieved when 
<fi 2 < V(cfi), that is, the potential energy of the scalar 
field dominates over the kinetic energy. Sustaining the 
accelerated expansion for long enough also requires the 
second derivative of the field (fi to be small. These con¬ 
ditions can be encoded in the smallness of two potential 
dependent slow-roll parameters: 


ev(< fi) = 



Vv(<fi) = M 2 Pl 



(49) 


Inflation proceeds when ey < 1 and \qv \ 1, and ends 

when ey ((fiend) ~ 1. During the expansion, the universe 
grows by a number of e-folds: 


JVM-log (5f )=/'"'m»£ 


1 


d (fi 


\j2ty Mpi 

( 5 °) 

The Cosmic Microwave Background fluctuations are cre¬ 
ated about 40 to 60 e-folds before the end of inflation [34] . 
The precise value depends on the details of reheating, 
post-inflationary thermal history and the energy scale of 
inflation. The integral constraint N((fi crn i,) « 40 — 60 
provides us with the field value when the CMB fluctua¬ 
tions are created. This in turn can be used to find the 
main observables related to the power spectrum of the 
fluctuations on the sky: 


A, 


V 


24 n 2 M^e v ’ 
n s — 1 ~ 2rjv — 6ey, 
r « 16ey, 


(51) 


where the scalar amplitude A s , the scalar spectral index 
n s , and the tensor-to-scalar ratio r are evaluated at (fi cm b- 

We see from Eqs. (49) and (51) that n s and r do not 
depend on the overall factor AT 4 . Scanning the interest¬ 
ing parameter ranges, —10~ 4 < E < 0 and 0 < c < oo, 
the primordial tilt n s changes by at most 0.5% and the 
tensor-to-scalar ratio r by a couple of per cent. Thus the 
slow-roll results stemming from the potentials considered 
here are fairly universal as well as completely consistent 
with the current observational constraints coming from 
Planck [35]. We plot the predictions for the ( n s ,r) pair 
in Fig. 2. 

On the plateau, the potential of the field is V ~ \E\ A/ 4 . 
As a result the amplitude of the scalar perturbations 
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FIG. 1. The dependence of the potential from Eq. (44) on 
the parameter c. Note that the shape does not change for 
c > 0(1). Here, we choose E ~ —10 -3 . 
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FIG. 2. Predictions for the ( n a ,r ) plane for — 1CF 1 < E < 0 
against the observations ( Planck TT+lowP data). 



A s ~ M 4 |A|/M 4 ; is degenerate in parameters E and 
M. Observationally, the amplitude is given by A s ss 
2.3 x 10^ 9 . For the ranges of E and c considered here, 
M can range from about two orders of magnitude below 
the GUT scale all the way to the Planck scale where our 
formalism breaks down as there is very little RG flow 
and gravitational effects become important. As a limit¬ 
ing case the potential in Eq. (44) gives the correct size 
of A s for the combination M ~ Mpi and E ~ —10 -1 °. 


C. The discrete case with E = — 1 


In the remainder of this paper, we will limit ourselves 
to the theory with E = — 1 in Eq. (38). As argued 
above, this special choice of the eigenvalue corresponds 
to the IR theory that is connected to the free theory. In 
terms of our quantum mechanical problem, this is the 
un-normalisable solution that contains the lowest-state 
normalisable mode with a discrete spectrum of integral 
eigenvalues. Eq. (40) now gives us the potential for (fi , 






















FIG. 3. The potential from Eq. (52) depending on the pa¬ 
rameter c. Note that the shape doesn’t change for c > 0( 1) 

and that for small c a kink appears in the potential. FIG. 4. Predictions for the (n s , r) plane against the observa¬ 

tions (Planck TT+lowP data) when E *■ — 1. 


which is 


M 

M 4 


1 

2 



2 


- log 



I C(C) 
M 4 ' 


(52) 


We expect the overall scale M to be low, as we are ex¬ 
amining the dominant solution of the ERG in the IR. 
As in the general case, fixing the constant C makes the 
potential only depend on the overall scale parameter M 
and the shape parameter c. 

Again, at the minimum of the potential, d^V (< <t>min ) = 
0 , the mass of the scalar field is independent of the shape 
parameter c, 


dlV(tmin) = M 4 /M 2 Pl . (53) 

Furthermore, when c > 0(1), the potential and its 
derivatives do not depend on c and the potential has 
a stable shape (see Fig. 3). The plateau of the po¬ 
tential and the quadratic behaviour near the minimum 
are also fairly independent of c. The parameter governs 
how sharp the transition is between the two regimes and 
where it happens. For smaller c the transition happens 
at higher (j>. At c = 0 we restore the familiar quadratic 
potential of the free theory. 

In the plateau region we use the asymptotic expansion 
of the error function, 



which is approximately valid for <f> > AMpi. As before, 
in Eq. (46), we recover the logarithmic behaviour remi¬ 
niscent of the radiative corrections to a flat potential. In 
the special case considered here (E = —1), the correc¬ 
tions can be turned off for c = 0, when the free Gaussian 
theory is recovered. 

The slow-roll predictions of our model are plotted in 
Fig. 4 against the observational constraints. The spectral 


index and the tensor-to-scalar ratio are independent of 
the overall scale M. For the given shape of the potential 
(i.e. fixed c) the values of n s and r are determined. Then 
M is determined such that the amplitude of fluctuations 
matches the observations. Hence, the inflationary model 
described by the potential in Eq. (52) would fall in the 
class of ‘One parameter’ models [6]. 

For the ‘theoretically preferred’ values of the shape pa¬ 
rameter, c ~ 0(1), we are in very good agreement with 
the observations. Parameter c measures the deformation 
of the potential away from the free theory (c = 0). For 
the small values (10 -8 < c < 10“ 14 ) the scalar to tensor 
ratio r can be reduced by several orders of magnitude 
while being consistent with the observations. This re¬ 
sults in the further lowering of the scale of the potential, 
M, to about three orders of magnitude below the GUT 
scale. This provides more room for the RG flow and adds 
to the validity of our analysis. 

For the cases of small c the potential we study de¬ 
velops a kink (see Fig. 3) that introduces non-negligible 
derivatives of the field <p. This makes the slow-roll predic¬ 
tion less accurate. The transient violation of the slow-roll 
conditions in the potential can give rise to observable fea¬ 
tures in the primordial power spectrum. This motivates 
further study with the exact integration of the inflation¬ 
ary background and perturbations together with the re¬ 
laxation of the mean-field approximation in our deriva¬ 
tion of the potential. 


IV. CONCLUSION 

In this paper, we studied non-perturbative scalar field 
potentials derived as an IR limit of the exact renormalisa¬ 
tion group equation within the mean-field approximation. 
We demonstrated that this approximation is valid when 
considering the potentials of an almost constant field. 
This setup is precisely what is required to start the ex¬ 
ponential expansion of space - inflation. The potentials 
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we derived were capable of supporting slow-roll inflation 
with the values of the spectral index and tensor-to-scalar 
ratio fully consistent with the recent observations. De¬ 
spite our very general treatment of the scalar field po¬ 
tentials, the slow-roll results are largely independent on 
the constants parametrising the potentials and thus re¬ 
sult in fairly universal predictions. For the special case 
where we required that the IR theory was continuously 
connected to a free Gaussian fixed point, we identified 
a parameter range that led to transient violation of the 
slow-roll conditions. Hence, this results in a possibly ob¬ 
servable feature in the primordial power spectrum, which 
should be further studied in the future. 
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Appendix A: Decomposition of the UV theory in 
the mean-field approximation for scalar theories 


In this appendix, we consider in greater detail the de¬ 
composition of UV scalar theories in terms of the ERG 
equation basis presented in Eq. (22) . We will assume the 
space-time derivatives of (j> to be very small, so that we 
can use the mean-field approximation in which we can 
treat 4> as a constant real variable x. In this approxima¬ 
tion, the path integral (19) becomes a one-dimensional 
integral 


Z= f V(j>9 uv [0]« f dxe~ s ^ x \ (Al) 

J J —OO 


where the Euclidean action S(x) equals the potential 
V(x) of the UV scalar theory. 

In order for the expression (Al) to be well-defined, we 
assume that the theory is normalisable, i.e. that the inte¬ 
gral over x in (Al) is finite, cf. Eq. (42). It is important 
to note that in a general quantum field theory, i.e. in 
the absence of a precise definition of the path integral, 
this need not be the case. All that is usually assumed is 
that the potential is bounded from below. As a simple 
example of a normalisable theory, we can think of the 


theory, which gives 




dx 


— km x — 4 t\x 





1/4 



(A2) 


In the expression, K a (x) is the modified Bessel function 
of the second kind. In the small A and m expansions, 


4^ +0 ( A2/ "‘ 9 )' < A3 > 

Z *‘ =r (i)(^) 1/4+ °( m2/A3/4 )' (M) 

respectively. The expansions imply that at m = 0, a 
perturbative treatment around A = 0 would be ill-defined 
as the partition function is divergent. 

As an example of a well-known and well-defined QFT, 
which is “un-normalisable” from the point of view of our 
analysis, one can think of the Sine-Gordon potential. In 
this case, 


Zsg 


dx exp 




(A5) 


is divergent and a more general approach would be 
needed to treat such theories in the MFA of the ERG. 
In such cases, the simple procedure of Eq. (Al) is insuf¬ 
ficient for computational purposes. 

To show how the computation of the ERG-basis de¬ 
composition (22) can be performed in simple example of 
scalar UV theories, we restrict our attention to normal¬ 
isable theories with polynomial potentials and consider a 
single scalar field theory in the MFA with 

S( x) = ^ m 2 x 2 + U(x). (A6) 

We assume that S{x) gives a finite Z, as for example in 
the </> 4 theory case computed in (A2), for which U{x) = 
Aad/4!. 

The general solution to the ERG equation (21) was 
presented in Eq. (41). For completeness, we restate it 
here: 


\H(a;) =e = x ~ 


Co iF\ 


(1 + E. 1 x 2 \ 

V 2 ; 2 ; T J 


+ D 0 x 



(2 + E_ 3 x 2 \ 
V 2 ’2’T) 


(A7) 


To form a complete set of functions out of (A7), we will 
restrict ourselves to the basis of normalisable solutions for 
which E are negative integers, i.e. E £ {—1, —2, —3,...}, 
with Go = 0 for E £ {-2,-4,...} and D 0 = 0 for 
E £ {—1, —3,...}. Since the first parameter of the hyper¬ 
geometric function is now always an integer, it is useful to 
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rewrite the confluent hypergeometric functions in terms 
of the Hermite polynomials: 


i-Fi 



(—l)"ra! 
(2n)! 


II> 



(A8) 


V2x\F\ 



(—l)"n! 
(2n + l)! 


H‘2n+l 



(A9) 

The Hermite polynomials multiplied by exp{—a; 2 /4}, i.e. 
the functions in the set of 


2 ' 4 H n 



n > 0, 


(A10) 


form a complete orthonormal basis on the Hilbert space 
of L 2 -integrable functions on the entire real axis x £ 
(—00,00). The completeness relation has the form' 


dxe~ x2 / 2 H n 




V2n2 n n\ S nm , 
(All) 


where n,m £ {0,1,2,...}. Using the Hermite polynomial 
basis, we can now write the ERG flows of normalisable 
T uv theories as 

Vuv{t,x) = . (A12) 


Note also that ^uv{x) = ^uv(t = 0, x). Using the com¬ 
pleteness relation (All), we can express the coefficients 
7 n of the 'teuv expansion (computed at t = 0) as 

1 " = v^!L dx<>uv{x)H "(.7i)- (A13) 


In our prototypical example of the </> 4 theory, for which 
'S'uv is even in x (and y), it is clear that only with 
even n contribute to f(/y. Similarly, this would be true 
in any theory with an even U(x). To solve for the 
coefficients in such theories, it is particularly convenient 
to write the Hermite polynomials in terms of associated 
Laguerre polynomials, 


H 2 n{y) = (—l)"2 2rl n! L,fj -1 / 2 ) ( y 2 ) , 
H 2n+1 (y) = (—1)"2 2n+1 n\L£M 

since they obey the recurrence relation 


i!."i (y 2 ) = 1 


^ 2n + 1 + a — y 2 


n+ 1 
n + a 


(y 2 


(A14) 

(A15) 


< a “) 


7 For a detailed discussion of Hermite and Laguerre polynomials, 
their orthogonality relation and the proof of completeness, see 
Chapter V. of the reference [36]. 


Note that y = x/y/2. Using Eq. (A16) with a = —1/2 
and a = 1/2, we can find the recurrence relations for 
all even and odd 7„ from the knowledge of (70,72) an d 
(71,73), respectively: 


1 


7n+2 = 


(n + 2 )(n + 1) 


dln_ 

dm 2 


( 

1 

(" + b 

7 n + ^7n-2 


(A17) 


To see how this works on a specific example, let us 
again return to the case of the </ 4 theory, for which 


7o = 


3m 2 3mi ( 3to 4 


1 


72 = 


37T 


8 V m 2 A 3 


V 4A )' 

(3TO 4 + (2 + TO 2 ) A) I1/4 


-TO 2 (3to 2 + A) /_!/4 


+3 TO 4 L 


'5/4 


3to 4 

"4A" 


- I 


3/4 


( 3TO 4 


V 4A 


(A18) 



(A19) 


The rest of the coefficients 7„ can then be generated using 
(A18). Note that l a (x) is the modified Bessel function 
of the first kind. As a check on the results, we note that 
in the UV (at t = 0), we need to integrate over the entire 
real axis, x £ (—00,00). In that case, only the n = 0 term 
contributes to the integral over the sum (A12). In the 
extreme IR limit of t —> 00, the dominant contribution 
to the RG flow of Tj/y(t,a;) comes from the n = 0 (or 
E = —1) term, i.e. the Gaussian effective action. 

We note that the effective IR potential of the <f> 4 theory 
is not of the form we used to study inflation. This is be¬ 
cause in this appendix we restricted ourselves to only the 
basis of hypergeometric functions for which the first pa¬ 
rameter was an integer. In case of the </ 4 theory and other 
normalisable theories, this is sufficient to decompose the 
entire theory and solve the RG flow. However, we expect 
that for more general theory, in particular those with 
un-normalisable potentials, this type of a decomposition 
would not suffice and a more general basis would be re¬ 
quired. 


Appendix B: Momentum Space Solution 

In this appendix, we solve the system in momentum 
space, where the solution takes a much simpler form. We 
first recall the equation 

dt'b(x) = Hty(x) = (dl + xd x ) \H(a;), (Bl) 

where t is the appropriately rescaled RG time. Note 
that the Hamiltonian is symmetric in x, so that so¬ 
lutions decompose into symmetric solutions ’I's}#) and 
anti-symmetric solutions 4/ a (a:), respectively. 
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Let us go to momentum space and write 
tf(s) = f *(p)e ipx , 

J P 


We can now expand ’L(g) as usual in a basis {i/jk = 

e~ Ekq }, 


(B2) 


§(<?) = f C(E) e~ Eq . (BIO) 

Je 


where 



Note that 'L(g) is closely related to the Laplace transform 
of C(E). Indeed, we will see below that we need C(E) = 
0 for E > 0. The basis {ipk = e~ Ekq } is orthogonal with 
respect to the inner product 


Plugging this into Eq. (Bl), we find 


C01,^2) 


i p'LOC 

— / d qe q ^~ E ^ = 5{Ei - E 2 ). (Bll) 

27TZ J —ioQ 



^(p) +P 2 ^{p) - P)pd P ) 


e ipx = 0. 


(B4) 


This gives the RG equation in momentum space after an 
integration by parts on the last term, 


Note that in terms of p , restricting ourselves to E £ Z, 
the expansion (BIO) is just the regular Laurent expansion 
of ’L(p). 

We can now solve the re-defined RG equation (B9). 
That is, we assume that the UV theory has the expansion 


9f^ip) = - {p 1 2 3 + 1 + P d p ) T (p). (B5) 


*uv(p) = f C{E)p~ E . (B12) 

J E 


Note that the equation is a first-order differential equa¬ 
tion. Let us perform a redefinition of coordinates, 

pd p = d q , (B6) 

which is satisfied by 

p = Ke q , (B7) 


The solution of (B9) then reads 

*t(p)= [ C{E)e tE p~ E . 
Je 

The “eigen-theories” are therefore given by 


'J'(p) = 


e 2y 


(B13) 


(B14) 


for some constant K. We can set K = 1 without loss of 
generality. The RG equation then reads 

d^{q) = -(e 2q + l + d q )*(q). (B8) 

Rescaling the theory as 

T(g) = exp f^e 2q + q) T(g), 


As t —>■ oo, we expect that the eigen-theories with largest 
E, for which C{E ) ^ 0, will dominate. Fourier trans¬ 
forming (B14) back to position space and taking the sym¬ 
metric part then gives 


^ s (a;)oce 1 F 1 


fl + E' 1 x 2 \ 

V 2 ; 2’T ) ’ 


while taking the anti-symmetric part gives 


T a ( x)(xxe \Fi 


2 + E. 3 a; 2 \ 

2 ; 2 ; y)’ 


(B15) 


(B16) 


we find 


d^{q) = ~dg^{q). 


as expected. Moreover, we find that the inverse transform 
is ill-defined for E > 0. That is, we require that the UV 
(B9) theory has C{E) = 0 for E > 0. 
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